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The spin-orbit splitting of the electron levels in a two-dimensional quantum dot in a perpendicular 
magnetic field is studied. It is shown that at the point of an accidental degeneracy of the two lowest 
levels above the ground state the Rashba spin-orbit coupling leads to a level anticrossing and to 
mixing of spin-up and spin-down states, whereas there is no mixing of these levels due to the 
Dresselhaus term. We calculate the relaxation and decoherence times of the three lowest levels due 
to phonons. We find that the spin relaxation rate as a function of a magnetic field exhibits a cusplike 
structure for Rashba but not for Dresselhaus spin-orbit interaction. 



I. INTRODUCTION 



I Recent years have seen an increasing interest in the 
spin properties of nanostructures^. Manipulation and 
readout of spins in solids could open the way to the de- 
velopment of a generation of electronic devices such as 
spin transistors, spin filters, and spin memory devices. 
In addition, the spin of an electron confined to a quan- 
tum dot (QD) is a promising candidate for a quantum 

[ bit^. Owing to the zero dimensionality of QDs, the 
electronic orbital states are quantized and the electron 
spin states are very stable due to a substantial suppres- 
sion of spin-flip mechanisms^. Progress in nanotech- 
nology has allowed the fabrication of QDs with desirable 
electronic and spin properties^sSiSifi. However, only re- 
cently it has been possible to measure the spin of an 

] electron in a QD. A single electron spin has been de- 
tected by magnetic resonance force microscopyiS and 
the readout of an individual electron spin in a QD via 
pulsed relaxation measurementsii and optical orienta- 
tion experiments^ have been reported. In these exper- 
iments, an external magnetic field was used to distin- 
guish spin-up and spin-down states split by the Zeeman 

I energy. Spin relaxation measurements between Zeeman 
levels in a QEiiLi^ confirm the theoretical predictions 
that spin-flip relaxation in a QD is suppressed with re- 
spect to a bulk structure^i^ii^. Indeed, very long single- 
spin relaxation times have been observed: up to 0.85 ms 
in two-dimensional (2D) GaAs QDs^-'^, and up to 20 ms 
in self-assembled GalnAs QD»i^. The spin relaxation 

[ is expected to be dominated by hyperfine interactions 
with the nuclei at magnetic fields below 0.5 Tii^iiSii^ 

' and by spin-orbit (SO) interactions for magnetic fields 
of about 1 Tesla (see Ref. |^ and for higher magnetic 
fields (see Ref. IT^ . In general, the SO interaction con- 
sists of two distinct contributions: the Dresselhaus SO 
coupling^^ which is due to bulk inversion asymmetry of 
the lattice and the Rashba SO couplingi^ which is due 
to structure inversion asymmetry along the growth di- 
rection. Both of these SO terms result in the splitting 
of electron energy levels and in the mixing of the elec- 
tron spin states. The latter makes spin-flip relaxation 



between Zeeman levels possible, for example, due to the 
phonon scattering. Note that usually it is not simple to 
separate these two SO mechanisms and estimate the rel- 
ative contributions of each SO term. In experiments, to 
obtain information about one of the SO couplings, nor- 
mally the other is neglectedi^iSSiSi. This leads to a lack 
of precision in estimates of the SO coupling strength and 
to a neglect of the effects of the interplay of the Rashba 
and the Dresselhaus SO couplingsi^iSSiS^. Hence, it is 
very important to find a way to separate these SO mech- 
anisms, to increase our understanding of the SO relax- 
ation processes, and to improve predictions of the spin 
properties of nanostructures. It is well known^ that for 
2D quantum wells the different SO couplings can be dis- 
tinguished experimentallji22i2S42& via detection of the as- 
sociated anisotropy of the spin splitting in the conduction 
band. In contrast, such a detection is not possible in QDs 
since the spin splitting of the levels, being quadratically 
in the SO coupling, is isotropic. Still, as we point out 
now, the SO couplings in QDs can be distinguished via 
their associated spin relaxation rates since they strongly 
differ due to different level mixing properties. 

In this paper, the electron energy spectrum and the 
spin relaxation for a 2D QD in magnetic fields perpendic- 
ular to the QD surface are studied. Level anticrossing— 
(due to the SO coupling), at a point of accidental level 
degeneracy (due to the interplay between the orbital and 
magnetic confinement), is analytically investigated. This 
anticrossing is caused by the Rashba SO term only, lead- 
ing to a cusp structure in the magnetic-field dependence 
of the spin relaxation rate, whereas the spin relaxation 
rate due to the Dresselhaus SO coupling is a monotonic 
function of magnetic field in this region. This qualita- 
tive difference in the spin relaxation for different SO cou- 
plings can serve to extract the different contributions in 
SO coupling. 



II. MODEL AND ENERGY SPECTRUM 

We consider a 2D isotropic QD with parabolic lat- 
eral confinement potential. An external magnetic field 
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is applied perpendicularly to the surface of the QD. The 
Hamiltonian of this system reads 

p2 1 1 

^0 = ^ + 2™*^o {^^ + y") + l^ai^sBG,, (1) 

where P = p+ (|e|/c)A(r), A(r) = (B/2)(-y, a;, 0) is the 
vector potential in the symmetric gauge, loq is the char- 
acteristic confinement frequency, and cr = {<Tx,(Jy,az) is 
the vector of the Pauli matrices. 

The SO interaction is taken into account by adding the 
linear DresselhausiiS^ and Rashbai^ terms for conduc- 
tion band electrons in a [001] two-dimensional electron 
gas (2DEG), 

Hd = /3(-a,P, + ayPy), Hr = a{a,Py - GyP,). (2) 

The axes x, y, and z are aligned along the principal crys- 
tallographic axes of GaAs. 

It is convenient to introduce new phase coordinates 
(qi, q2, pi, P2) which connected to the previous ones 
(x, y, px, Py) by the following formula^: 

P2 
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where 



n = 



Here lJc — \e\B/m*c is the cyclotron frequency. In the 
new phase coordinates, has the canonical form 



p\ 



pI 



"^' 2 2, 2 2^ 1 



2m* 



— K^i +^^2'72) + i^9l^BBaz. (3) 



The first level is the ground state. In the case of weak 
magnetic confinement (wq ^ (^c), the second level is 

lower than the third one (-Eool ^ ^iot)- However, at 
high magnetic fields, when the magnetic confinement is 
much stronger than the lateral confinement (wq <SC Wc), 

^ooT ^o/^c- The condi- 



^ool > -^S' because E^^^ 



tion for a crossing of the levels e'^I and SJqt given by 



r(0) 



OOi 



^lOT 



uji — ujz- In other words, this level crossing takes place 
when the magnetic length Ib = y/h/mFujc is equal to 

h[9*{g* + 1)]"^^^: where Iq = yJhJrnFojQ is the character- 
istic lateral size of a QD and g* = \g\m* /2mi^. Note that 
the level crossing occurs at accessible magnetic fields for 
QDs with lateral size > 15 nm. 

Now we take SO coupling into account and find the 
energy spectrum and eigenstates of electrons in a QD. 
For a GaAs QD the SO lengths are Ad = h/m*(3, Xr = 
h/m*a ~ 8 fiir?^ and are much larger than the hybrid 
orbital length I — ^h/m*Q, of a QD (Ad, \r ^ I). 
Therefore, the SO terms can be considered as small per- 
turbations. 

First we consider the Dresselhaus SO coupling [see 
Eq. It is important to note that in first-order 

perturbation theory there is no SO interaction between 
the levels E^^^ and E'f^^^ due to the Dresselhaus term 
((00 J, |i?D|10 t) — 0). Hence we can apply standard 
perturbation theory for nondegenerate levels. Thus, in 



first-order perturbation theory, we get En = En 



(0) 



|i) = |00T) + 



-|io i), 



\oo i) ~ ^^l^\oi T), 



UJ2 — UJZ 



^ iiot) + ^^!(!Z^l20i). 

oji + UJZ 



(4) 
(5) 
(6) 



Now we consider the Rashba SO coupling term. In 
this case, there is a SO interaction between the levels 
-EqqI and E^o'y Therefore, applying perturbation theory 
for degenerate levels, we have 



In this case, Hq can be considered as the Hamiltonian 
of two independent harmonic oscillators with hybrid fre- 
quencies Therefore, the energy spectrum and eigen- 
states of electrons in a QD without the SO coupling are 
given by 

E^°X^^ = huji{n + l/2) + hi02{m+l/2)-huJzs,, 
(qiq2\nmsz) = ^>„((7l\/"^*t^l/?l)$m(g2\/"l*W2/?l) |s^) , 



El ^ hfl- ^hojz, E2,3 = hn + ^{ujiT ujr) , (7) 



|1> - |00T>-A|0U>, 



where 



(8) 



|2) = cos||00i)-sin||10T)+Asin^|lU), (9) 
|3> = sin^|00i)+cos||10t>-Acos^|lU>,(10) 



where n,m = 0, 1, 2, . . ., Sz — ±1/2 is the electron-spin 
projection on the z-axis, uJz — IglfJ-BB/h is the Zeeman 
frequency, and <&«((?) are oscillator functions. 
Let us consider the three lowest levels: 



E, 



(0) 



^ hn- hujz/2, E^^l] = nn + hujz/2 



ooT 

_e|qi = hfi + huji — hujz /2. 



UJR 

tan 7 
A 



^{uji-uJzY+A{l/\Ryujl 
~2{1/\r)uji/{uji - UJz), 

[1/\r)uJ2/{uJ2 +UJz). 



(11) 



As can be seen from Eq. ((JJ), in the case of strong lat- 
eral confinement [uji — ujz ^ {1/Xr)uji\, E2 = E^^^ 
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and = e[^p but in the case of strong magnetic 
confinement [ujz — uji ^ (1/ X^) lui], the levels E2 and 
i?3 change places: E2 — i?|oT' ^3 ~ ^ool- 
crossing point for the levels i?ooi ^^'^ ^loj ("^i ~ '-^z)^ 
E'2.3 — + fiLUz/'2 =F ^r)^ z ■ Therefore, the Rashba 
SO coupling leads to an anticrossing of the levels E2 
and E3 at the point of accidental degeneracy of the lev- 



els E, 



(0) 



d E'iol [see inset in Fig. [Tl;b)]30. The dis- 
tance between the levels E2 and E^ at the anticrossing is 
A = 2{l/XR)hujz. For a GaAs QD with hcoo = 1.1 meV 
and Xr = 8 /zm, this anticrossing is too small for ex- 
perimental observation (A = 0.5 /xeV), but for an InAs 
(^g w li^ and Xr ~ 0.1 ^m-^) QD with the same size, 
the anticrossing can reach 0.1 meV. Note that this anti- 
crossing features were numerically studied for narrow-gap 
QDs in Ref.H^ 

Let us study the states |2) and |3). As can be seen from 
Eqs. © and ifuji^ujz > {l/Xn^i [7 = 0(1/ Xp)], 

|2) = |GG i) + 0{1/Xr), |3) = |10 T) + 0{1/Xr). 

With increasing B, the Zeeman energy becomes larger 
than hioi. In the case of loz — wi ^ {1/Xr)uji, 7 — 
TT + 0(1/ Xr) and these states change place. Therefore, 
the spin flips with a transition trough the anticrossing 
region. In the region of the anticrossing (7 w — 7r/2), 
the SO coupling of these states due to the Rashba term 
becomes essential and leads to a mixing of spin-up and 
spin-down states: 



|2) 
|3) 



(|00 D-l- 
(-100 i> 



|10 ^))/V2 + 0{1/Xr.), 
H10T»/V2 + OG/Ah) 



Note that, although for a GaAs QD the level anticross- 
ing is a quite small effect, the mixing of spin-up and 
spin-down states occurs in a sufficiently large region of 
magnetic fields (for a GaAs QD with Hloq ~ 1.1 meV and 
Xr = 8 /im, the mixing occurs essentially in the region 
of width « I Tesla) and thus can be observed exper- 
imentally. Indeed, let us consider relaxation processes 
between the state |1) (spin- up) and the states |2), |3). 
Beyond the mixing region one of the latter states is spin- 
up, the other is spin-down. Spin-flip relaxation is much 
slower than orbital relaxatior^, therefore, relaxation to 
the ground state from the states |2) and |3) is very differ- 
ent. However, in the region of mixing of spin-up and spin- 
down states, the spin-flip relaxation strongly increases 
and becomes comparable with orbital relaxation. Note 
that these anticrossing features in semiconductor QDs are 
very similar to the "hot spots" in polyvalent metalsS^. 

Moreover, it is interesting to note that spin relaxation 
due to the Rashba SO coupling differs from that due to 
the Dresselhaus SO coupling in this mixing region. As 
mentioned above, in the case of the Dresselhaus term 
there is no SO interaction between the states |2) and 
|3) [see inset in Fig. H^a)], therefore, there is no spin 
mixing of these states. Thus, spin relaxation due to the 
Dresselhaus SO coupling does not undergo a considerable 



increase, in contrast to spin relaxation due to the Rashba 
SO coupling. Note that, in the general case, when the 
SO coupling includes both the Rashba and Dresselhaus 
terms, there is no interplay between the Dresselhaus and 
Rashba terms in the spin relaxation rate in perpendicular 
magnetic fieldsi^ and the total rate is just the sum of two 
terms caused by these SO couplings. Therefore, we can 
study these two terms separately. 



III. SPIN RELAXATION 

We consider next phonon-induced relaxation in a QD. 
The coupling between electrons and phonons with mode 
kj (k is the phonon wave vector and j is the branch index 
j ~ L, Tl, T2 for one longitudinal and two transverse 
modes) is given hyi^ 



^/2pVks~/h 



(e^kj - ikE]^j)e 



ikr L + 



(12) 

where p is the crystal mass density, V is the volume of 
the QD, Sj is the sound velocity, A\^j = S,i£,id^ Pum, ^ = 
k/A:, d'^-' is the phonon polarization vector, Skj is the 
deformation potential, and Pum is the piezotensor, which 
has nonzero components only when all three indices i,l,m 
are different: P^yz = Pxzy = ■ ■ ■ = hn/es (es is the static 
dielectric constant). For GaAs, e/ii4 = 1.2 x 10'' eV/cm, 
es = 13.2, Skj = 5jxSo, and So = 6.7 eV). In Eq. ^ 
we introduced the form- factor Fikz) which is determined 
by the spread of the electron wave function in the z- 
direction: F{kz) = / dze''^''\'4)[){z)\'^ , where V'o(^) is the 
ground state envelope wave function of an electron along 
the z-direction. The form factor F{kz) equals unity for 
\kz \ ^ d~^ and vanishes for lA:^! » d~^ (see Ref. Il3l). 

Let us find contributions to the relaxation rate of tran- 
sitions between the levels |1) and |2) (r2i); |2) and |3) 
(r32); |1) and |3) (Fai). In the framework of the Bloch 
- Redfield theory, the phonon-induced relaxation rate 
(1/Ti) of a two- level system is a sum of transition prob- 
abilities between levels accompanied by absorption and 
emission of phonons^i and, for a QD, the decoherence 
time is T2 — 2T^. Therefore, using Fermi's golden rule 
and the expressions for the three lowest levels with the 
Dresselhaus SO couphng [Eqs. © and 0], we get the 
rates 



r2i 



l^Lul {N^^ + 1/2) f uji 



UJ2 



SnhpXjj + UJZ UJ2 — 0JZ 



^32 — 



W^P (Ny, + 1/2) / UJl 



UJ2 



32TrhpXjj \ll>i + u>z UJ2 ~ LOz 



(13) 



(14) 
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FIG. 1: (color online). Contributions to the relaxation rate l/Ti of phonon-induced transitions between the states |1), |2), 
and |3) of a GaAs QD with hiuo = 1.1 meV and d = 5 nm due to (a) the Dresselhaus and (b) the Rashba SO couplings 
(Ad = \r = 8 /im). The dashed and dot-dashed curves are orbital relaxation rates, the solid and dotted curves are the 
relaxation rates with a spin flip. The crossing (a) and the anticrossing (b) of the levels Ez and E2 are shown in the insets. 
The cusplike structure of the spin relaxation curve due to the Rashba SO coupling is caused by the mixing of the spin-up and 
spin-down states at the anticrossing. 



where w = uji — ujz, = ie^'^"/^ — l) ^ , and 



27r pTi/2 

"'0 



16) 



In the case of parabohc confinement along the growth 
direction of a QD, I^"^\x) can be expressed in terms of 
error functions (see APPENDIX EJ. 

In the case of Rashba SO coupling alone, we have 



(jv^_ + 1/2) f . {1/\r)oJ2 ,„ 
r2i = ^ 77^ ^ sm7/2H cos 7/2 



E 



%-Kpm*Q, \ " UJ2+UJZ 

H-/2s] A3) 



X > e 



(17) 



327rp(TO*S2)^ \2 UJ2+UJZ 



j 

{N^^ + 1/2) f . 
lpm*n ^ (^^»^7/2-sm7/2 



(18) 



{l/\li)uj2 
UJ2 + UJz 



(19) 



where w± = (wi + ujz i ljji)/2 and is defined by 
Eq. mi)- 



IV. ANALYSIS AND DISCUSSION 

Figure ^ shows these contributions to the relaxation 
rate due to the Dresselhaus and the Rashba SO couplings. 
As can be seen from this figure, the orbital relaxation rate 
(the dashed and dot-dashed curves) is independent of the 
SO coupling. The behavior of the spin relaxation rate 
(the dotted curves) is qualitatively the same for both the 
Dresselhaus and the Rashba SO couplings. Solid curves 
correspond to the the spin relaxation rate T^'^ between 
the Zeeman-split orbital ground state levels (F^^^ = F21 
for the Dresselhaus SO coupling and in the case of the 
Rashba SO coupling F^^^ — F21 on the left side of the cusp 
and T^'^ = F32 on the right of the cusp). Significantly, 
F^1^, in the case of the Rashba SO coupling, possesses a 
cusplike structure at the anticrossing point?—, whereas, in 
the case of Dresselhaus SO coupling, T^^ is a monotonic 
function of B^^. 

It should be noted that at B > 1 T the relaxation 
due to deformational acoustic (DA) coupling is much 
faster than that due to piezoelectric (PE) coupling, ex- 
cept in the case of orbital relaxation at high magnetic 
fields, when relaxation induced by PE-phonons is of the 
same order as that due to DA-phonons. Since d <ti Iq 
and q sa l^^, the factor F{qz) « 1 in the electron- 
phonon interaction operator (see Ref. Il3j) and the re- 
laxation is practically independent of d aside from the 
orbital relaxation at low magnetic fields: V3i{B = 0) « 
^o(^<^o + l/2)Sge-'^o<iV2s? /67rpm*sJ (the spin relaxation 
rates are zero at i? = 0). The orbital relaxation rate has 
a maximum when the phonon wave length is comparable 
to the lateral size I of a QD {ql « 3). At high mag- 
netic fields, the orbital relaxation rate decreases with 
B [as {luq/lOcY' for DA coupling and as {loq/lOcY for 
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PE coupling], since uji — > Wq/wc at high B. The rate 
oc oj'^(uji — ujz)^ at low magnetic fields, is zero at 
the anticrossing, and oc at high magnetic fields. 
The spin relaxation rate between the Zeeman-split levels 
(X ijj^ (at low magnetic fields k — 7 for DA coupling 
and fc = 5 for PE coupling. At high B, k = i for DA 
coupling and k — \ for PE coupling). In the anticrossing 
region, the spin relaxation rate due to the Dresselhaus SO 
coupling is a monotonic function of B: V^^ = oc lo\, 
but that due to the Rashba SO coupling has a strong 
increase at the anticrossing point and near this point 
F^T cx w|/[(l-wzM)^+4(Z/Afl)2]. Therefore, there is 
both a qualitative difference (in the magnetic-field depen- 
dence) and quantitative difference (at 4.8 T the Rashba 
SO coupling gives V'^'^ « 10^ s~^ but the Dresselhaus SO 
coupling gives V^'^ ~ 70 s^^) in the behavior of the spin 
relaxation rate F-'-^ due to the Dresselhaus and Rashba 
SO coupling. This can serve as a means of extracting in- 
formation on the different contributions to the total SO 
coupling strengths^. 

Note that, with a decrease in the lateral size of a 
QD, the cusp and the maximum in the orbital relaxation 
rate are shifted to high magnetic fields. For a larger SO 
coupling (smaller SO length), the spin relaxation rates 
have higher values, because F^T ^-^2^ ^j^g (.^gp 
shape is smoother. The temperature dependence of the 
relaxation rates is only important for transitions between 
the levels with a separation comparable to the tempera- 
ture: the rates decrease with temperature for the orbital 
relaxation at high magnetic fields (when the level spacing 
~ hw'^/ujc), for the spin relaxation between Zeeman-split 



levels (when the Zeeman energy ~ T) at low magnetic 
fields, and for F32 at the anticrossing (when the level 
spacing \uji ~ ujz \ ^ T). 



V. CONCLUSIONS 

We have shown that at an accidental degeneracy point 
the Rashba SO coupling leads to an anticrossing. The 
mixing of the spin-up and spin-down states at the an- 
ticrossing enhances the spin relaxation rate due to the 
Rashba SO coupling relative to the spin relaxation rate 
due to the Dresselhaus SO coupling. 
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APPENDIX A: PARABOLIC CONFINEMENT 
ALONG THE ^-DIRECTION 

In the case of parabolic confinement along the growth 
direction of a QD, F{kz) = exp{~d?k'^/4:), where d is the 
width of the quantum well, and integrals in Eq. (|15|l can 
be expressed in terms of the imaginary error functions 
erfi(a;). After some algebra we get 
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where t| = {P - d'^) /2s^j. 
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